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_i ■ Abstract. This paper presents a study of the accretion properties of 19 very low mass objects (M*~ 0.01 — 0.1 Mq) 

in the regions Chamaeleon I and p Oph. For 8 objects we obtained high resolution Ha profiles and determined 
mass accretion rate M ac and accretion luminosity L ac . Pa/3 is detected in emission in 7 of the 10 p Oph objects, 
but only in one in Cha I. Using objects for which we have both a determination of L ac from Ha and a Pa/3 
detection, we show that the correlation between the Pa/3 luminosity and luminosity L ac , found by Muzerolle et 
al. ( 1998) for T Tauri stars in Taurus, extends to objects with mass ~ 0.03 Mq; L(Pa/3) can be used to measure 

^ ■ Lac also in the substellar regime. The results were less conclusive for Br7, which was detected only in 2 objects, 

\£) ' neither of which had an Ha estimate of M ac - Using the relation between L(Pa/3) and L ac we determined the 

accretion rate for all the objects in our sample (including those with no Ha spectrum), more than doubling 
the number of substellar objects with known M ac - When plotted as a function of the mass of the central object 
together with data from the literature, our results confirm the trend of lower M ac for lower M*, although with a 
large spread. Some of the spread is probably due to an age effect; our very young objects in p Oph have on average 
an accretion rate at least one order of magnitude higher than objects of similar mass in older regions. As a side 
product, we found that the width of Ha measured at 10% peak intensity is not only a qualitative indicator of the 
accreting nature of very low mass objects, but can be used to obtain a quantitative, although not very accurate, 

C^h' estimate of M ac over a large mass range, from T Tauri stars to brown dwarfs. Finally, we found that some of our 

>■» k objects show evidence of mass-loss in their optical spectra. 

Cu • 1. Introduction In this paper we discuss one important parameter in 

£> this scenario, namely the accretion rate in very low mass 

Several objects of very low mass discovered in regions of objects, which, together with the disk mass, determines 

star formation have infrared exc ess typ ical of circu mstel- the timescale for disk survival. Mass accretion rates have 

^ lar disks (e .g., N atta and Testi |200IJ Testi et al. |2002J been determined for a handful of very low mass objects, 

Natta et al. !2UU2p and sho w sign s of accretion-relat ed ac- mostly from model-fitting of the Ha profiles (Muzerolle 

tivity (Jayawardana et al. 12002| White and Basri |2003| et al. I2UU51. The results indicate that the trend of lower 
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Muzerolle et al. [2003): they undergo an evolutionary accretion rates for objects of lower mass, known for TTS 

phase similar to t hat of classical T Tauri stars (TTS) (White and Ghez I27M1 Rebull et al. ESTTfTllj) continues in 

(Jayawardana etal.HSSSl- This evidence sets potentially the substellar regime, where typical accretion rates are 

important constraints on the formation of very low mass mucn i ower than in TTS and can be as low as - 5 x 10~ 12 

objects, and on the dominance of dynamic al pro cesses in M Q /yr (Muzerolle et al. l2UB0) l. There is also evidence that, 

star formation in general (e.g., Bate et al. [20031). which, aga i n as f or TTS (e.g. Calvet et al. l2UB0) l. the fraction of 

however, need to be made more quantitative. Statistics accreting objects among very low mass objects is higher in 

of disks, measurements of their masses, characterization younger regions (Jayawardana et al. l2H0ll20031jjl . All this 

of the accretion as a function of the mass of the central po ints to a continuity in the accretion properties across 

objects in regions of different age and properties are re- the range of mass explored so far, which suggests a similar 

quired, and are being accumulated, exploiting fully the formation mechanism from solar to sub-stellar masses. 

capabilities of 8-m class telescopes. r™ i, u u j c 

Inese results, nowever, are based on lew measurements 

only, and need to be confirmed. In particular, it is impor- 

Send offprint requests to: natta@arcetri.astro.it tant to extend them to younger, more embedded regions, 

* Based on observations collected at the European Southern where the use of Ha to derive the accretion rate is made 

Observatory, Chile. hard or altogether impossible by the high extinction. In 
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these regions, it is natural to use instead infrared hydro- 
gen recombination lines, such as Pa/3 or Br7. For TTS, 
it has been shown that the flux in these lines provides a 
good estimate of the accretion rate (Muzerolle et al. 11998(1 . 
The purpose of this paper is to extend this method to the 
sub-stellar regime, by obtaining near-IR spectroscopy of a 
sample of well-known objects in two star-forming regions 
(Chamaeleon I and p Oph) for which we have also acquired 
and analyzed high-resolution Ha profiles. 

The range of masses of interest goes from the low end of 
the stellar range (<0.2M Q ) to brown dwarfs (BD: <0.075 
Mq) and planetary-mass objects, below the deuterium- 
burning limit (<0.013 M© or 13 Mj). In fact, these distinc- 
tions are meaningful only "a posteriori" , when bona-fide 
stars reach the main sequence, not in young star-forming 
regions, where all objects, independently of their mass, 
derive their energy from contraction. For the sake of sim- 
plification (and in view of the uncertainties affecting mass 
determinations), we will talk in the following of very low 
mass objects or VLMOs. 

In §2 we detail the observations, obtained with UVES 
and ISAAC on VLT. The results are presented in §3 and 
discussed in §4. A summary follows in §5. 

2. Observations 

2.1. Sample 

We have collected a total of 19 objects, 9 in Chamaeleon I 
and 10 in p Oph. Their properties are listed in Table 1. For 
the Cha I sample, spectral types, effective temperatures, 
extinctions, luminosities and masses are from Comeron et 
al 1(2000(1 . Quoted uncertainties are ±150 K in T e ff, ±0.2 
dex in Log L*, and a factor of 2 in M*. Stellar parameters 
for the p Oph sources are from Natta et al. (2002), with 
the exception of GY10, not included in the Natta et al. 
sample, for which we adopt the parameters of Wilking et 
al. {1999). Uncertainties are ±100 K in T eff , ±20-30% in 
L*, ± 20-40% in M* (see Testi et al. fMffl for a detailed 
discussion). The sample contains a number of bona-fide 
BDs, several objects that are either very low mass stars or 
relatively high-mass BDs, one object (p Oph-033) which 
has an estimated mass of ~10 Mj only. Based on their 
location on the HR diagram (see Fig. QJ, the VLMOs in 
Cha I are significantly older than those in p Oph (>2 Myr 
versus <1 Myr, respectively). 

2.2. Low resolution near-infrared spectroscopy 

Near-infrared low resolution spectra in the J and K bands 
were obtained for all our targets using the ISAAC near- 
infrared camera and spectrograph at the ESO-VLT UT1 
telescope. The observations were carried out partly in 
Visitor Mode on April 8, 2002, and partly in Service Mode 
from March to July 2003 (see Table 2; Table 2 is only avail- 
able in electronic form at "http://www.edpsciences.org"). 
We used the 0.6 arcsec slit and the low resolution grisms 
that offered a ~800 resolution across the wavelength 
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Fig. 1. Location of the sample objects in the HR diagram. 
Filled dots are objects in p Oph, open squares objects in 
Cha I. The tracks are from D'Antona and Mazzitelli 1(1997(1 
and their unpublished 1998 update. 



ranges 1.0-1.3 pm and 1.95-2.45 /im. The observing se- 
quence was composed by a set of pairs of spectra with 
the target in different positions along the slit, to allow for 
an efficient and accurate sky subtraction. The on-source 
integration times varied from 25 to 40 minutes depend- 
ing on target brightness and observing band (the J-band 
observations were typically 40-30% shorter than the K- 
band ones). Standard calibrations (flats and lamps) and 
telluric standards spectra were obtained for each obser- 
vation. Wavelength calibration was performed using the 
lamps observations and adjusted by measuring a few OH 
sky lines before sky removal. We did not attempt to obtain 
flux calibrations. 



2.3. High resolution optical spectroscopy 

High resolution optical spectra for the eight brightest (in 
the optical bands) objects in our sample were collected 
using the UVES spectrograph on the ESO-VLT UT2 tele- 
scope (see Table 2). For each target a series of three to 
five ~45 min exposures were obtained in the period from 
March to June 2003. All spectra were obtained using the 
standard UVES setup with central wavelength 580nm. 
The chosen setup allowed us to cover a nominal spec- 
tral range from ^450 to 680nm. However, due to the very 
low signal obtained in the blue part of the spectra (450- 



Natta et al.: BD accretion 



580 nm), only the red part was analyzed. The slit used 
was 12 arcsec long and 1 arcsec wide, giving a resolu- 
tion R~40000. Standard calibrations were obtained by the 
ESO staff as part of the Service Mode and the data were 
reduced using the UVES pipeline. All spectra were of ex- 
cellent quality and the Ha line was detected with high 
signal to noise in all individual spectra. 

3. Results 

3.1. Accretion rates from Ha 

All eight observed objects have Ha emission, in seven 
cases with equivalent width >10 A. Of these, five (2 in 
Cha I and 3 in p Oph) show in all our observations broad 
Ha profiles. This, as discussed by White and Basri (2003) 
and Jayawardana et al. (2002J), is indication of accretion- 
related rather than chromospheric activity. For each of 
these objects, we have computed an average profile, shown 
in Fig. [3 and compared it to the prediction of the mag- 
netospheric accretion models developed by Muzerolle et 
al. (2001) and recently applied to a sample of sub-stellar 
objects by Muzerolle et al. (I2U03|I . Accretion rates derived 
from Ha profiles agree well with those obtained from veil- 
ing in TTS (e.g., Muzerolle et al. HM5|l and in the few 
VLMOs for which M ac could be estimated from veiling 
measurements (Muzerolle et al. [2003 ) . 

The underlying assumption is that the entire Ha (as 
well as the other hydrogen lines) comes from accreting 
matter, and that the contamination from wind/jet emis- 
sion is negligible. This is not necessarily true in all strongly 
accreting TTS (see, for example, Bacciotti et al. I2()()2|l . 
but is probably a good approximation for our objects (see 
§3.4). 

For all objects we have adopted a stellar mass M*=50 
Mj, radius 0.5 R© and a magnetospheric truncation ra- 
dius of 2.2-3 R*. The results are not very sensitive to the 
exact values of these parameters. The fitting procedure 
gives for each object a value of M ac and of the inclination 
angle i. The resulting accretion rates are given in Tabled 
they range from 10~ 9 M Q /yr for p Oph-102 to - lO" 11 
M Q /yr for p Oph-030. The inclination angle is ~ 65 — 75 
deg for Cha Ha2, p Oph-030, p Oph-032, - 80 deg for Cha 
Ha6,and ~30 deg for p Oph-102. The fits are of acceptable 
quality, with the exception of Cha Ha2, whose broad and 
flat profile (clearly present in all the observations) cannot 
be reproduced by the current models (similar profiles are 
observed in few other VLMOs and discussed in Muzerolle 
et al. l2003fl . The fit can be somewhat improved by adopt- 
ing an ad-hoc temperature profile for the accretion flow, 
whereby the temperature has been reduced relative to the 
fiducial model in the outer part of the flow near the disk. 
This results in less emission near the line center, the most 
uncertain part of the profiles given the approximations 
used in the radiative transfer calculations (see Muzerolle 
et al. l2001[) . Nevertheless, both models for Cha Ha2 have 
similar inclination and an accretion rate of the order of 
10~ 10 M Q /yr, which we will adopt in the following. As 
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Fig. 2. Normalized average Ha profiles for the 5 objects 
with broad lines. The solid curves are the observed pro- 
files, the dashed ones are the model fits. For Cha Ha2, 
the dotted lines shows a model with an ad-hoc temper- 
ature profiles (see text). The stellar continuum has been 
subtracted. 



discussed by Muzerolle et al. (2003), the values of M ac 
derived from the Ha profiles should be accurate within a 
factor of ~ 3 — 5. 

A second group of 3 objects, all in Cha I, (Cha Hal, 
Cha Ha3, Cha Ha5) has always narrow Ha profiles, typi- 
cal of cromospheric activity. For these objects, we estimate 
an upper limit to M ac of about 10 -12 Mo/yr, for which 
the magnetospheric accretion models predict a line emis- 
sion below the detection limit. 

In Fig. we compare our estimates of M ac to those 
of other VLMOs and TTS from the literature. The figure 
plots Mac as a function of the full width of Ha measured at 
10% of the peak intensity (after continuum subtraction). 
Dots are VLMOs (mass <0.2 M Q ), where M ac is measured 
by fitting the Ha profiles (filled dots from this paper, open 
dots from Muzerolle et al. 2003); pentagons are VLMOs 
with values of M ac high enough to produce measurable 
veiling in the optical (White and Basri 2003; Muzerolle 
et al. I2TJUU1 I2HU51 Barrado y Navascues et al. l2lMjl . The 
squares show the location on this diagram of a sample 
of more massive TTS, from which M ac has been derived 
from veiling (Gullbring al. I1998[) . and the 10% Ha width 
are from high resolution Ha profiles obtained by Edwards 
et al. (1994). The total sample spans a range of masses 
from about 0.04 to 0.8 M . Fig. |3| shows that our deter- 
minations of M ac agree with the general trend shown by 
other objects from the literature, and that there is conti- 
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Table 1. Object Properties and Hydrogen Lines 



Star 


ST 


T e ff 


LogL* 


Av 


M* 


Hq10% 


EW(Hq)° 


EW(Pa/?) a 


EW(Br 7 ) a 


LogM a 6 c 


LogM a c c 






(K) 


(Lq) 


(mag) 


(Mj) 


(km/s) 


(A) 


(A) 


(A) 


(M /y) 


(M Q /y) 


Cha Hal 


M7.5 


2770 


-1.96 


0.2 


40 


155 


35 


<0.3 


<i 


< -12 


< -11.4 


Cha Ha2 


M6.5 


2910 


-1.47 


0.8 


70 


366 


33 


0.3 


<1.5 


-10.0 


-10.8 


Cha Ha3 


M7 


2840 


-1.46 


0.3 


60 


142 


10 


<0.9 


<1.2 


< -12 


< -10.1 


Cha Ha5 


M6 


2980 


-1.31 


1.0 


100 


137 


6.5 


<0.9 


<1.7 


< -12 


< -10.0 


Cha Hq6 


M7 


2840 


-1.57 


0.26 


50 


274 


48 


<0.3 


<1 


-10.5 


< -10.8 


Cha Hq7 


M8 


2690 


-2.19 


0.3 


30 


- 


- 


<0.6 


<1.0 


- 


< -11.1 


Cha Ha9 


M6 


2980 


-2.26 


1.5 


70 


- 


- 


<0.6 


<1.0 


- 


< -11.4 


Cha HalO 


M7.5 


2770 


-2.28 


0.1 


40 


- 


- 


<0.3 


<0.9 


- 


< -12.0 


Cha Hall 


M8 


2690 


-2.44 


0.0 


30 


- 


- 


<0.3 


<0.9 


- 


< -12.1 


p Oph-023 


M7 


2650 


-1.40 


8.0 


40 


- 


- 


1.80 


<1 


- 


-9.3 


p Oph-030 


M6 


2700 


-1.15 


3.0 


60 


292 


30 


0.30 


<1 


-10.8 


-10.1 


p Oph-032 


M7.5 


2600 


-1.22 


2.0 


40 


248 


50 


0.40 


<0.9 


-10.5 


-9.8 


p Oph-033 


M8.5 


2400 


-2.10 


7.0 


10 


- 


- 


<0.7 


<2.1 


- 


< -10.3 


p Oph-102 


M6 


2700 


-1.10 


3.0 


60 


377 


40 


2.0 


<1 


-9.0 


-8.9 


p Oph-160 


M6 


2700 


-1.40 


6.0 


45 


- 


- 


3.3 


4.0 


- 


-9.0 


p Oph-164 
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6.0 
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0.80 


<1.6 
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-9.3 


p Oph-176 
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- 
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<0.5 


<1.3 
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< -9.7 


p Oph-193 


M6 


2650 


-1.00 
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60 


- 
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1.8 


2.9 


- 


-8.9 


p Oph-GYIO 


M8.5 


2350 


-1.44 


18 


30 


- 


- 


<0.5 


<1.0 


- 


< -9.0 



a Positive values indicate emission. 
6 Accretion rates from Ha profiles. 
c Accretion rates from Pa/3 luminosities. 



nuity between measurements obtained from veiling (which 
are fully independent of the Ha width) and those derived 
from Ha profiles. 

In the sub-stellar regime, the Hal0% width, which 
can be derived directly from the observations, is consid- 
ered a good indicator of accretion, with the separation 
between accretors and non-accretors set at Hal0% width 
between ^200 km/s ( Jayawardana et al. 2003b) and ~ 270 
km/s (White and Basri 2003). Fig. [21 confirms this result. 
Adopting ^200 km/s as a limit, one would missclassify one 
accreting object out of 23 (MHO-4, with M oc ~ 1.5 x 10" 11 
M /yr and Ha 10% width of 154 km/s; Muzerolle et 
al. HJnH and two non accretors (V927 Tau and USco CTI- 
075) out of the 37 for which there are M ac upper limits 
(Muzerolle et al. 2003 and this paper). Furthermore, Fig. [3] 
shows that there is a rather good correlation between the 
Ha 10% width and M ac over the whole range of mass 
from BDs to TTS, so that it is possible to use the ob- 
served width not only to discriminate between accretors 
and non-accretors but also to get an approximate estimate 
of the accretion rate, without performing detailed model 
fits. We show in Fig.|3]the best-fit relation between these 
two quantities for Ha 10% width > 200 km/s, which can 
be expressed as: 



LogAf ac = -12.89(±0.3) + 9.7(±0.7) x 10~ 3 Hal0% (1) 



where Hal0% is the Ha 10% width in km/s and M ac is 
in M Q /yr. 

The spread is rather large, and can be due in part 
to the fact that in most cases the measurements used to 
derive M ac and the high-resolution Ha profiles have not 
been obtained simultaneously. We show, as an illustration 
of possible problems, the rather unusual case of the TTS 
DF Tau which has an accretion rate of about 10 -7 M Q /yr, 
based on veiling from medium-resolution spectrophoto- 
metric data obtained in 1996 (Gullbring et al. I1998fl . 
High resolution profiles of several Balmer lines obtained 
non-simultaneously from 1988 to 1990 by Edwards et 
al. (1994), show broad emission in all the lines, but with 
large discrepancies between them; in particular, the Ha 
width is smaller than that of the higher Balmer lines. One 
could certainly improve the correlation if larger and simul- 
taneous sets of data for TTS were available. M ac values 
derived in this way are necessarily inaccurate for individ- 
ual objects, and should be used with care. Nevertheless, 
they can be very useful when dealing with large samples 
of objects. 

3.2. Paf3 and Brj as accretion indicators 

Our near-IR spectroscopy does not allow us to resolve the 
line profiles, and we can only measure equivalent widths, 
listed in Table 1. Of the 9 Chamaeleon I objects, Pa/3 



Natta et al.: BD accretion 





-6 








i 


1 1 1 | 


I 1 I 1 I 

/ 
/ _ 

■ 

r 

DF Tau 

■ 1"" 

■ -" 


o 
S 

o 
o 
to 

s 


-8 












- 


ao 


-10 










o o o« 


- 















9 ^° 


. 










O 




o • 
















o 


_ 




-12 




OOD 


taomm 


DO 
1 


o 

I 


lii" 



200 400 600 

Ha 10% Width (km/s) 



Fig. 3. Accretion rate as a function of the Ha full width 
at 10% peak intensity. Dots (full and open) and pentagons 
are objects with mass <0.23 M@. Full dots are from this 
paper, open dots from Muzerolle et al. 1)2003): for both 
samples M ac has been derived from model-fits of the Ha 
profiles. Pentagons are substellar objects for which M ac 
has been obtained from veiling measurements (Muzerolle 
et al. 120011 120001 White and Basri 120051 Barrado y 
Navascues et al. 2004). The small dots at Log M ac =— 12 
are VLMOs with no evidence of accretion (37 in total); 
the value of the accretion rate is the model-determined 
upper limit. Squares are TTS (mass > 0.3 M©) for which 
accretion rates have been measured from veiling. Most are 
TTS in Taurus- Auriga (M ac from Gullbring et al. 110051 
Ha widths from Edwards et al. 1994). We have also added 
TW Hya (Muzerolle et al. l2000|l and GN Tau from White 
and Basri (2003). The dashed line shows the best-fitting 
relation (Eq. 1). The four connected squares show the 10% 
full width of Ha H/3, H 7 and H<5 for DF Tau from Edwards 
et al. (|1994[1 : Ha is the narrowest of the four. 



is detected in only 1. For the others, we can set upper 
limits to the equivalent width (3<r) of about 0.3-0.9 A. In 
Ophiucus, on the contrary we detect Pa/3 emission in 7 of 
the 10 observed objects. Br 7 is detected only in 2 objects 
in p Oph. Of the 5 objects with evidence of accretion from 
Ha, Pa/3 is detected in 4, while Br 7 is not detected in any. 
Fig-Elshows the spectra of the objects with line detections. 
These results suggest that Br 7 , which is intrinsically 
weaker than Pa/3, is more difficult to detect even in highly 
reddened objects, where differential extinction favors lines 
at longer wavelengths. For this reason, we will focus in 
the following on the use of Pa/3 as an accretion indica- 
tor. We have therefore converted equivalent widths into 
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Fig. 4. Near-IR spectra of objects with Pa/3 (top panel) 
or Br 7 (bottom panel) detection. 

fluxes using broad-band magnitudes and then correcting 
for extinction for the Cha I objects, while for the p Oph 
objects we have used the calibrated low-resolution spectra 
of Natta et al. (|2002[1 for the continuum and corrected for 
extinction the resulting line fluxes. 

Fig. shows the relation between the Pa/3 luminos- 
ity (derived from the measured fluxes assuming a distance 
of 150 pc for Ophiucus and 160 pc for Cha I) and the 
accretion luminosity, computed from M ac and the stel- 
lar parameters given in Table 1. According to Muzerolle 
et al. I|1998f) . these two quantities (rather than line flux 
versus M ac ) show the tightest correlation in TTS, and 
we have added Muzerolle's sample to ours. The figure 
shows that the trend of lower Pa/3 luminosity for lower 
L ac extends to the very low values of L ac which character- 
ize VLMOs. In fact, one can derive a rather good relation 
to estimate L ac from the Pa/3 luminosity over the whole 
range of masses from few tens of Jupiter masses to about 
one solar mass, shown by the dashed line in Fig. [3] and 
given by: 



LogL ac = 1.36(±0.2) x LogL(Pa /3 ) + 4.00(±0.2) 



(2) 



where L ac and L(Pa/3) are in units of L Q . Note that this re- 
lation is, within the uncertainties, identical to that derived 
by Muzerolle et al. (£1998) for TTS over a much narrower 
interval of L(Pa/3). 
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Fig. 5. Pa/3 luminosity as a function of L ac . The dots 
are the VLMOs from this paper. The squares are T Tauri 
stars in Taurus- Auriga from Muzerolle et al. Ijl998|l . The 
dashed line is the best-fitting relation (Eq. 2). 
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As for Ha we have ignored the possibility that some 
Pa/3 emission may come from wind or jets, driven by ac- 
cretion. In several TTS, the profiles of Pa/3 and Br7 are 
well described by magnetospheric accretion models (Folha 
and Emerson I200ip . but there are some examples of Pa/3 
wind emission (Whelan et al. 2004). Inspection of Fig. 
shows that, very likely, any correction for a wind/jet con- 
tribution to the Pa/3 luminosity will be lost within the un- 
certainties of the individual measurements and the large 
scatter of the points. Nevertheless, the possibility of a con- 
tribution from outflowing matter to lines should be kept 
in mind when discussing individual objects. 

Using the relation (2), we have derived L ac and M ac 
for all the objects in our sample; M ac is given in the last 
column of Table 1. The largest differences between this 
determination of M ac and that from the Ha profile is of 
the order ±0.7-0.8 in log, i.e., a factor ~6 (see also Fig.|3J). 

3.3. Variability 

Spectroscopic and photometric variability over a large in- 
terval of timescales is typical of pre-main sequence stars of 
all mass, including sub-stellar objects. By necessity, this 
introduces an artificial spread in correlation such as those 
of Fig. [3] and [5] when the accretion rate or luminosity and 
the line properties are not measured simultaneously. We 
expect that the quality of these correlations could improve 
significantly if simultaneous observations were available. 



Another consequence of variability is the difficulty of 
assigning a well-defined accretion rate to any specific star. 
Variations of the Br7 flux are typically within a factor 2- 
3 (see, for example, Greene and Lada 119571 Luhman and 
Ricke l*H?9"5)l : if similar variations apply to Pa/3 (for which 
we could not find multiple-epoch data in the literature) 
this would change the estimate of L ac by a factor ~ 2 — 5. 
In some cases, the lines become undetectable, and one 
would then classify the star among non-accretors. 

Similar uncertainties are obtained if we consider the 
variations of the Ha 10% width. In our spectra (see Fig.[J)J), 
the largest variation is seen in p Oph-102, where the Ha 
10% width varies from 335 km/s on May 29, 2003 to 411 
km/s on June 4, 2003. According to eq.(l), the correspond- 
ing accretion rate would increase by a factor ~5. In all 
the other objects the Ha variations are smaller, implying 
a maximum change of M ac of a factor ~3 at most. A com- 
parison of the sub-stellar objects in IC 348 for which Ha 
profiles have been obtained by Muzerolle et al. ( 2003 ) and 
Jayawardana et al. I|2003b|l shows variations that would 
results in changes of L QC again by a factor ^2 — 5, accord- 
ing to Eq.(l). 

An interesting case is that of p Oph-030 (GY 5), which 
has a strong Ha emission in our and Jayawardana et 
al. H2002f) spectra (but a very different profile and a 10% 
Ha width of 292 and 352 km/s respectively) but is nar- 
row (177 km/s) and weak in the spectrum of Muzerolle et 

al. <msi). 

Some of these variations may be due to differences in 
the quality of the spectra (resolution, signal-to-noise etc.), 
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since all the VLMOs are very weak, often close to the in- 
strumental limits. Still, most of the variability is likely 
to be real, as suggested by very similar results in TTS. 
For any individual object, one should have high quality 
spectra and long time sequences, which could provide not 
an instantaneous but a "typical" accretion rate. However, 
if we consider a large sample of objects, these individual 
variations will average up and will only increase the spread 
of any existing correlation, but not change it systemati- 
cally. 



3.4. Other Optical Lines 

The UVES spectral range contains a number of other 
emission lines (in addition to Ha), usually considered good 
indicators of mass-loss and accretion. Table [3] gives the 
equivalent width of the lines (positive values for emission) 
and, in some cases, as for the NaD doublet, some infor- 
mation on their profile. The spectrum of Cha Hal is too 
noisy, and we have dropped this object from the table. 

The Hel lines are generally interpreted as emitted in 
the accretion shock at the base of the infalling matter: 
these lines are indeed present and measurable in all the 
objects where we detect accretion in Ha (Cha Ha2, Cha 
Ha6, p Oph-030, p Oph-032 and p Oph-102), but not in 
Cha Ha3 and 5, which have narrow Ha profiles and no 
evidence of accretion. This result is consistent with our 
Ha analysis. 

Some of the objects show indications of outflows. These 
are particularly strong in p Oph-102, which has emission 
in all the forbidden lines of [01], [Nil], [SII], and shows a 
P-Cygni profile in the NaD doublet (see Fig.EJ), but also 
p Oph-32, Cha Ha2 and Cha Ha3 have detectable emis- 
sion in some of the forbidden lines, and clear emission in 
the NaD doublet. All these objects are good candidates 
for further studies of winds and jets in VLMOs. However, 
there is no evidence of H 2 vibrationally excited emission in 
the IR spectra of our objects, with the possible exception 
of p Oph-033, where, however, it is likely due to contami- 
nation from VLA 1623 (see also Williams et al. [1995). 

Weak accretion-driven winds are expected in VLMOs, 
if they behave as TTS. Note, however, that even in p Oph- 
102 the equivalent width of Ha is much larger than that of 
the [SII] lines (the ratio is ~ 200). The weakness of the for- 
bidden lines in comparison to Ha in our objects supports 
the assumption that the hydrogen lines form mostly in the 
accreting gas (see also Fernandez and Comeron 2001 and 
references therein). 

As well known the presence of Li is a strong indication 
of substellar nature and/or youth. The covered spectral 
range includes the Li I doublet at 6708 A. We detect it in 
absorption in all but one objects with equivalent widths 
ranging between 0.2 and 0.5 A, confirming, in combination 
with the late-M spectral types, that they are members of 
the clouds. The Li line is not detected in Cha Hal; we 
believe this is due to the very low S/N of the spectrum 
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Fig. 7. NaD lines in p Oph-102. 

which shows virtually no continuum (see also Jorgens and 
Giinther l2TM)l . 

4. Discussion 

4.1. Accretion as a function of mass and age 

Our sample of VLMOs more than double the number of 
objects with measured M ac . We plot in Fig.|S|our results, 
together with a compilation of data from the literature, 
as a function of the mass of the central object. The data 
include TTS and VLMOs, spanning a range of masses 
of about two orders of magnitude. All the masses have 
been redetermined self-consistently using the same evolu- 
tionary tracks (1998 unpublished update of D Antona and 
Mazzitelli 1997), and the accretion rates have been scaled 
to these values, when necessary. Note that we have plotted 
upper limits for M ac only for objects from this paper. 

Both M ac and M* have large uncertainties, discussed 
in detail, among others, by Muzerolle et al. (|2003|l and in 
this paper. There are, however, some clear trends. First of 
all, Fig. |H1 confirms, over a large range of masses, the trend 
already found by other authors (Muzerolle et al. [2003 
Rcbull et al. I200 0J) of decreasing M ac for decreasing M*. 
If, for example, we compare the median M ac for M*<0.1 
Mq to that in the mass interval 1-0.3 M Q , we find values 
of ~ 3 x 10~ 10 and ~ 10 -8 M /yr, respectively (neglect- 
ing upper limits to M ac ). It is unlikely that this trend is 
caused by sensitivity limits. Although, as we will discuss in 
the following section, in many VLMOs the accretion rate 
is close to the model sensitivity, this is not the case of 
TTS (see, e.g., Muzerolle et al. l2"UU3"j) . so that the correla- 
tion we find neglecting upper limits is probably somewhat 
shallower than the true one. 
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Table 3. Equivalent width of other optical emission lines 



Object [OI]A6300 [OI]A6363 [NII]A6583 [SII]A6716 [SIIJA6731 HeIA5875 HeIA6678 NaIA5896,5890 a 

(A) (A) (A) (A) (A) (A) (A) 

Cha Ha2 = - 

Cha Ha3 
Cha Ha5 
Cha Hq6 

p Oph-030 
p Oph-032 
p Oph-102 2: 0.4: 



0.06: 


0.06: 


<0.1 


1.6 


0.6 


E(PCyg) 


0.03: 


<0.1 


0.08 


<0.1 


<0.1 


E 


<0.1 


<0.08 


<0.1 


<0.2 


< 0.08 


E 


0.07 


<0.08 


<0.08 


1.1 


0.3 


E 


<0.1 


< 0.08 


<0.08 


1.0 


0.3: 


E 


0.1: 


0.12: 


<0.1 


6.0 


2.2 


E 


0.1: 


0.07 


0.2: 


4.4 


2.3 


E(PCyg) 



a E: both components in emission; PCyg: additional blue-shifted absorption (see Fig. |7J . 



A second effect is clearly shown in Fig. El namely that 
the accretion rate, for the same mass of the central object, 
is higher in younger regions. The p Oph brown dwarfs have 
accretion rates of roughly 10 -9 -10 -10 M Q /yr, more than 
an order of magnitude larger than most objects of similar 
mass in Taurus, Chamaeleon and IC 348. An age depen- 
dence of Mac is known in TTS (e.g., Calvet et al. l2TJ0"tjjl . 
and is suggested for VLMOs by the fact that the fraction 
of objects with narrow Ha profiles (non-accretors) is larger 
in older star forming regions (Jayawardana et al. 2002 



20033). Our results confirm this trend: only 10% (1 out 
of 10) of the objects in our Cha 1 sample are accreting; if 
we include also the Gomez and Persi (|2001[1 VLMOs, then 
the fraction becomes 28% (5 out of 18 objects). 



The dependence of M ac on age makes it hard to deter- 
mine the exact relationship of M ac vs. M* from samples 
where objects in regions of very different age are collected 
together, as in Fig.|5] For example, the p Oph BD popula- 
tion selected by Natta et al. (20021) IS particularly young, 
and it would be interesting to compare their accretion 
rates to equally young TTS in p Oph, rather than to TTS 
in Taurus. Also, there are a few VLMOs in Taurus, Cha 1 
and IC348 with very large accretion rates, comparable or 
even higher than those of the p Oph VLMOs. Some of 
these estimates need to be confirmed, and the possibility 
of having detected a flare ruled out. Still, one wonders if 
they could be much younger than the average population 
of the star-forming region to which they belong. 



In spite of the large uncertainties that affect individ- 
ual objects, these preliminary results show that it is now 
possible, with the advent of 8-m class telescopes that give 
access to the VLMOs population, to quantify the depen- 
dence of M ac on mass and age, overcoming, due to the 
large range of M* that one can explore, the uncertainties 
of individual measurements. For this one needs larger and 
more homogeneous samples of stars (not only VLMOs but 
also TTS) in a variety of star-forming regions. 
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Fig. 8. Accretion rate as a function of the mass of the 
central object. Filled dots and squares are VLMOs in 
p Oph and Cha I, respectively, from this paper; M ac has 
been determined from L(Pa/3) as described in §3.3. Arrows 
are 3<7 upper limits. Open squares are objects in Cha I 
from Gomez and Persi (2001) with Pa/3 emission, where 
M ac has been determined from L(Pa/3). Objects in IC 
348 (Muzerolle et al. 20031) are shown by diamonds, in 
TW Hya (Muzerolle et al. I2UU0|I by pentagons, in Taurus 
(White and Ghez ESUoTl White and Basri|j003, Muzerolle 
et al. I2003[) by crosses. We do not plot objects from the 
literature with no detected accretion. 



4.2. Disks and Accretion 

Classical TTS show a clear correlation between accretion 
and the presence of a circumstellar disk. Disks are de- 
tected around a number of VLMOs, which show mid- 
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Table 4. Correlation between accretion and mid-IR ex- 
cess 



Star 


Ace. 


Ace. 


Ace. 


6.7 pm 


14.3 /im 




(Ha) 


(Pa/3) 


(Br 7 ) 






Cha Hal 


N 


N 


N 


Y 


Y 


Cha Ha2 


Y 


Y 


N 


Y 


Y 


Cha Ha3 


N 


N 


N 


Y 


N 


Cha Hq5 


N 


N 


N 


Y 


N 


Cha Hq6 


Y 


N 


N 


Y 


N 


Cha Ha7 


- 


N 


N 


N 


N 


Cha Hq9 


- 


N 


N 


Y 


Y 


Cha HalO 


- 


N 


N 


N 


N 


Cha Hall 


- 


N 


N 


N 


N 


p Oph-023 


- 


Y 


N 


Y 


Y 


p Oph-030 


Y 


Y 


N 


Y 


Y 


p Oph-032 


Y 


Y 


N 


Y 


Y 


p Oph-033 


- 


N 


N 


Y 


Y 


p Oph-102 


Y 


Y 


N 


Y 


Y 


p Oph-160 


- 


Y 


Y 


Y 


Y 


p Oph-164 


- 


Y 


N 


Y 


Y 


p Oph-176 


- 


N 


N 


Y 


Y 


p Oph-193 


- 


Y 


Y 


Y 


Y 


p Oph-GYIO 


- 


N 


N 


Y 


Y 



infrared fluxes well in excess of the photospheric ones 
(e.g., Comeron et al. IT3M1 Persi et al. 12710711 Bontcmps 
et al. 1200 l[l . modeled as the emission of disks heated by 
radiation from the central objects (Natta and Testi I200"T1 
Testi et al. 12001 Natta et al. [20021 • 

Table 0] summarizes the information on the presence 
of a mid-IR excess in our sample objects. All 10 objects 
in p Oph have been detected by ISOCAM at 6.7 and 14.3 
pm (Bontemps et al. I2001JI . and the correlation between 
disk and accretion, i.e., broad Ha and/or emission in the 
hydrogen IR lines, is basically confirmed by the results of 
this paper. 

The Chamaeleon sample is in a way more intriguing. 
Three objects (Cha Hal, Cha Ha2, Cha Ha9) have excess 
emission at 6.7 and 14.3 /im(Persi et al. 2000), modeled 
by Natta and Testi (J2D01) as due to disks. Of these, Cha 
Ha2 is clearly accreting, Cha Hal and Cha Ha9 (for which 
we have no Ha profile) do not. Three objects have been 
detected at 6.7 /im, but not at 14.3 /im; of these, Cha 
Ha3 and Cha Ha5 do not show evidence of accretion in 
Ha nor in Pa/3 or Br7, while Cha Ha6 has broad Ha, 
but no emission in the IR lines. Finally, three (Ha7, HalO 
and Hall) have no detected mid-IR excess and no IR line 
emission. 

It appears from these results that disks are detected 
more easily than accretion, and that there is a popula- 
tion of VLMOs with evidence of disks but no detectable 
accretion activity. This is somewhat different from TTS, 
where most WTTS have no evidence of accretion but also 



no detectable IR excess (as some of the Cha I objects). 
This is probably due to the fact that even in very active 
VLMOs the accretion rate is close to the sensitivity limit 
of the various methods one can use, making it difficult 
to discriminate between "accretors" and "non-accretors" . 
On the contrary, in TTS typical accretion rates are at 
least two orders of magnitude higher than the sensitiv- 
ity limit of veiling measurements and/or Ha profile fitting 
(Gullbring et al. IT555} Muzerolle et al. IMtSj) . and the 
separation between accreting and non-accreting objects is 
therefore much more meaningful. 

One of the implications of the very low accretion rates 
that characterize BDs is that their disks should live long. 
If we assume, for example, that the disk mass scales with 
the mass of the central objects in BDs as in TTS, we have 
typical disk masses of about 10 -3 M Q ; with an accretion 
rate of ~ 3 x 10~ 10 M Q /yr, as measured in p Oph, the 
disk will live at least 3 x 10 6 years, as a typical TTS. This 
agrees with the results of Liu et al. (2003), who found no 
signficant difference in the fraction of disks between BDs 
and TTS, and that the age of the TTS with disks was 
similar to that of TTS without disks. 

It is more difficult to use this same result to rule 
out "truncated" disks, as those produced by the dynam- 
ical ejection of stellar embryos predicted, e.g., by Bate 
et al. 12003 ■ Let us assume, for example, that the disk 
has an outer radius of 100 AU, and that, once ejected, 
loses all material outside a radius of about 10 AU. The 
mass left depends critically on the unknown density pro- 
file of the disk, and may range from about 1/3 (if the 
surface density profile was S oc i?~ 15 ) to about 1/10 (for 
S ex R~ 10 ) of the original mass. A disk with mass of only 
10~ 4 M Q will live about 3 x 10 5 years, much less than a 
typical disk of TTSs, and also much less than the age 
of many BDs with infrared excess in star-forming regions 
(e.g., Jayawardhana et al. 2003a). However, if the density 
profile is very peaked, the difference will not be significant, 
given the large uncertainties and spread of all the param- 
eters involved. Only very accurate model predictions can 
be verified or falsified by this kind of arguments. 

A final point worth mentioning is that the relation of 
the accretion rate with the mass and age of the central 
object may shed light on the physical process by which 
disk matter accretes onto the central star. The simplest 
suggestion is that the accretion rate may be related to the 
temperatures in the disk, which in turn are closely coupled 
to the luminosity of the central object. This may explain 
both the dependence of M ac on M* and its decrease with 
time, since objects of all masses, and BDs in particular, 
fade as they get older. However, it is not clear that the 
change in temperature is big enough to account for the 
changes in M ac . More subtle effects can be at work, related 
for example to the level of chromospheric and coronal ac- 
tivity (and hence X-ray fluxes) of the star and resulting 
disk ionization (see discussion in Muzerolle et al. |2I 
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5. Summary and conclusions 

We present in this paper new optical and near-infrared 
spectroscopic observations of a sample of 19 very low mass 
objects in the regions Chamaeleon I and p Oph. Our main 
aim was to test if accretion rates could be determined 
from the intensity of the hydrogen recombination lines in 
the infrared, namely Pa/3 and Bi'7. This possibility is very 
important when dealing with large numbers of VLMOs 
in young star-forming regions such as p Oph, where high- 
resolution spectroscopy in the visual is possible for only a 
handful of objects. 

To this purpose, we have obtained high-resolution Ha 
profiles of all the objects in our sample bright enough to 
be detectable with UVES on VLT. The comparison of 
the observed profiles with those predicted by magneto- 
spheric accretion models (Muzerolle et al. 2001) shows 
that five objects are actively accreting, with M ac in 
the range ~ 10 -9 — 10 -11 M /yr, while the remaining 
three in Cha I have narrow and symmetric profiles, and 
Mac<W- 12 M /yr. 

Pa/3 is detected in emission in 7 of the 10 p Oph ob- 
jects, but only in one in Cha I. We show that the corre- 
lation between the Pa/3 luminosity and the accretion lu- 
minosity L ac , found by Muzerolle et al. (1998) for TTS 
in Taurus over a range of masses ~0.3-l M , extends to 
masses about ten times lower. Whereas the relation be- 
tween L(Pa/3) and L QC (eq.(2)) can certainly be improved 
when more data will be available, it is already reason- 
ably well defined, and can be used to measure L ac , and 
consequently M ac , for large samples of obscured VLMOs. 
The results are so far less conclusive for Br7 for which we 
have a detection in only 2 objects, none with measured 
accretion rate. 

Using eq.(2) and the measured L(Pa/3), we determined 
Mac for all the VLMOs in our sample, increasing the num- 
ber of VLMOs with known M ac by more than a factor of 
two. When plotted as a function of the mass of the cen- 
tral object together with all the existing data on VLMOs 
and TTS from the literature, our results confirm the trend 
of lower M ac for lower M+, although with a large spread. 
Some of the spread is likely due to an "age" effect, namely 
that, for a given value of M*, M ac decreases strongly with 
time. Our very young VLMOs in p Oph have on average 
M ac at least one order of magnitude higher than objects 
in older star forming regions. The dependence of M ac on 
central mass and time needs to be better constrained by 
future observations. However, it is already clear that these 
data can provide very valuable information on the ac- 
cretion mechanism in disks during the pre-main-sequence 
evolutionary phase. 

One side product of our analysis is that the width 
of Ha measured at 10% peak intensity correlates quan- 
titatively with M ac over a large range of M ac and M*. 
This supports the currently used criterion to separate ac- 
cretors from non-accretors among VLMOs according to 
the Ha 10% width (White and Basri ?M)^ and confirms 
that a value of 200 km/s (as suggested by Jayawardana et 



al. 12003 bfl is a rather good threshold value. Furthermore, 
it suggests that one can use the Ha 10% width to ob- 
tain a quantitative estimate of M ac over a large range of 
masses. Although the spread of the correlation (eq.(l)) is 
large, and there are a number of caveats to keep in mind, 
the possibility of estimating M ac without going through 
rather time consuming model fits is interesting, and de- 
serves further attention. 

Finally, we would like to note that high-resolution 
near-infrared spectroscopy on 8-m class telescopes is be- 
coming available. As with Ha for visible objects, high 
quality infrared line profiles can be compared to model 
predictions to derive accretion rates for obscured ones. 
The additional advantage of the IR lines is that they are 
likely less contaminated than Ha by wind/jet emission, 
flares and stellar activity in general; a sufficiently large 
sample of objects with measured Ha and Pa/3 profiles will 
be very valuable in providing more "calibration" points 
for the line intensity-accretion rate relations, and in bet- 
ter constraining magnetospheric accretion models. 
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